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We present a time-dependent Langevin description of dynamics of stock prices. Based on a simple sliding- 
window algorithm, the fluctuation of stock prices is discussed in the view of a time-dependent linear restoring 
force which is the linear approximation of the drift parameter in Langevin equation estimated from the finan- 
cial time series. By choosing suitable weighted factor for the linear approximation, the relation between the 
dynamical effect of restoring force and the autocorrelation of the financial time series is deduced. We especially 
analyze the daily log-returns of S&P 500 index from 1950 to 1999. The significance of the restoring force to- 
wards the prices evolution are investigated from its two coefficients, slope coefficient and equilibrium position. 
The new simple form of the restoring force obtained both from statistical and theoretical analyses suggests that 
the Langevin approach can effectively present the macroscopical and the detail properties of the price evolution. 

PACS numbers: 89.65.Gh, 05.45.Tp, 87.23.Ge 



I. INTRODUCTION 

The analysis of financial data by methods developed for 
physical systems has a long tradition |T, and has at- 

tracted the interest of physicists. One of the most motivated 
reasons is that it is a great scientific challenge to understand 
the dynamics of a strongly fluctuating complex system with a 
large number of interacting elements. In addition, it may be 
possible that the experience gained in the process of study- 
ing complex physical systems might yield new results in eco- 
nomics. 

There are many observables generated from financial mar- 
kets, and one central issue of the research on the dynamics 
of financial markets is the statistics of price changes which 
determine losses and gains. The price changes of a time se- 
ries of quotations x{t) are commonly measured by returns: 
r :— x{t + t) / x{t), log-returns, or increments: Ax := 
x{t + r) — x{t) at a time scale t. In 1900, Bachelier pro- 
posed the first model for the stochastic process of returns — an 
uncorrected random walk with independent identically dis- 
tributed Gaussian random variables However, prices do 
not follow a signal random walk process SQP, 8|. For ex- 
ample, the daily correlation has been known as the daily log- 
returns correlated with themselves in such a way that positive 
returns are followed by positive returns as well |9, 10|. Many 
considerations have been aroused by this effect recently and 
meanwhile the related research has been reported, not only 
for daily data fTHl but also for high frequency data fT2lfT3ll . 

The Langevin equation (LE) which distinguishes the de- 
velopment of sample path into the deterministic and random 
terms, has been used to deal with the Brownian motion prob- 
lem. Recently, the Langevin approach was used to analyze 
the financial time series on scale [ 14. 151 . Friedlich et al. 
lll4ll have investigated how price changes Ax on different time 
scales T are correlated motivated by hierarchical structure of 
financial time series, which is similar to the energy cascade 
in hydrodynamic turbulence. They derived a multiplicative 
Langevin equation from a Fokker-Planck equation (FPE) in 



the variable scale r and performed the statistical way to dis- 
tinguish and quantify the deterministic and the random influ- 
ence on the hierarchical structure of the financial time series 
in terms of the drift and diffusion parameters, D'^' and D^'^\ 
respectively 1 14]. 

Different from the former study in which the LE is used 
to analyze the scales evolutions of finance, yet, in this paper, 
with the Langevin description, a new insight on the dynam- 
ics of the process will be obtained by investigating the time- 
dependence of log-returns. The time-dependent properties of 
prices evolution are derived in the way of estimating drift pa- 
rameter A{z) of sampled local periods in the sliding window. 
Then, the relation between A{z) and autocorrelation C, aver- 
age return (z), from which the practical significance of A{z) 
can be recognized, are resulted both from the statistical time- 
dependence of A{z) and some theoretical analyses. Besides, 
our Langevin description contains, as a particular case with 
flat-A(z), the effect of daily correlation in log-returns. On the 
other hand, the form of diffusion parameter B{z) got in this 
paper, to some extent, explains the heavy tailed probability 
densities of price changes. 

The research are mainly carried out from the samples of the 
daily log-returns z{t)\T=id of S&P 500 index from 1950 to 
1999, containing 12583 days, thus covering a wide time range 
with many different economic and political situations. 

The paper is organized as follows. In Sec. II, taking the 
daily log-returns as an example, we generally discuss the ap- 
plication of Langevin approach to log-returns series. In Sec. 
Ill, we show the results and discussions. Finally, the summary 
and the outlook of this paper are given in Sec. IV. 



II. THE LANGEVIN APPROACH TO LOG-RETURNS 
SERIES 

For a time series of prices or market index values x{t), the 
log-return z{t) = z{t)r over a time scale r is defined as the 
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forward change in the logarithm of x{t), 

z{t)r = lnx{t + T)-lnx{t). (1) 

The behavior of daily log-return z (i) as a stochastic variable 
is described by the following LE: 

dz = A{z)dt + B{z)dw, (2) 

where the drift parameter A{z) and diffusion parameter B{z) 
respectively describe the deterministic and the random influ- 
ences on the time process of log-returns, and dw denotes the 
increment of a standard Wiener process. It is assumed that, 
within each sampling window the parameters may depend on 
the log-returns, but not explicitly on time (stationary). Thus, 
the drift and diffusion parameters of the sampled period can 
be extracted from the sampled data by simply using the defi- 
nition 1 16|, 

B{zf = ^im i^((Z(t + At)-z)2)U(,)^,. (4) 

Here (• • ■) denotes the averaging operator and Z{t + At) is 
a realization of the LE From Eq.(|3j, it is obvious that 
the drift parameter A{z) is the average increment of unit time 
under the condition Z{t) — z, which represents the determin- 
istic influences. B{z) is the deviation of A{z) which pictures 
the random influences. It has been known that the autocor- 
relation of the log-returns decays very fast which is usually 
characterized by a correlation time much shorter than a trad- 
ing day 1 5]. When the time increment A< is larger than 1 day, 
the daily log-returns can be considered as the result of many 
uncorrelated 'shocks'. Thus, in this paper. At is mainly set as 
1 day. Compared to the length of time window T, At = Id 
approximately accords with the limit in Eq.0 and Eq.@. 

Based on the samplings all over the long series, the statis- 
tical results of A{z) and B{z)'^ which are respectively esti- 
mated by Eq. ^ and Eq. have their simple and general 
forms. The results of A{z) are close to a linear form, and that 
of B{z)^ are close to a parabolic form, 

A{z) = az + b (5) 
B{zf = a'z^ + b'z + c'. (6) 

Fig. [l]presents the statistical results of A{z) and B{zY (cir- 
cles) which are estimated from the daily log-returns time se- 
ries of S&P 500 index from 28 May. 1962 to 10 Jun. 1970 
(window length T = 2000(i). The statistical results for large 
z are more noisy and uncertain than the points near the ori- 
gin, because these border points are visited rarely by the tra- 
jectory. On the contrary, as viewed from statistics, the more 
one given Zi is visited, the more times the averaging oper- 
ator (• • ■)\z(t)=zi in Eq. (|3jl works, which would produce 
more accurate and reasonable A{zi). Thus, while approxi- 
mating A{zi) (i = 1,2, ...,n) with linear form and B{zi)'^ 
(i = 1, 2, n) with parabolic form by Least-squares fit, the 
effect of the visited probability for each Zi should be consid- 
ered, with each Zi corresponding to its own weight. 
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FIG. 1: (Color online) The statistical results of A(z) and B{z) (cir- 
cles) estimated from the daily log-returns of S&P 500 index from 
28 May. 1962 to 10 Jun. 1970, a window with length T = 2000 
days. The approximations of the statistical results with weight fac- 
tor P{zi) (solid lines) and that with equal-weight (dashed lines) re- 
spectively correspond to the fitting errors a a = 2.564 X 10~* and 
8.650 X 10"^ and crs = 1.336 x 10"^ and 5.144 x 10"*. The 
expressions of solid lines are, A{z) = -0.7569a - 3.24305 x 10"* 
and B{z) = 0.92478^^ - 0.005622 + 2.48829 x 10"^ 

It is natural for a physical scientist to define the weighted 
factor of Zi as its probability, 

P(z,) = Niz,)/T, (7) 

where the N{zi) is the frequency of Zi within the time window 
of length T. Fig.^shows the approximations of the statistical 
results of A{z) and B{z)'^ with equal-weight (dashed lines) 
and weighted factor P{zi) (solid lines). The values of the fit- 
ting error <ja and as (the mean standard deviations from ap- 
proximations of A[z) and B{z)^) with weighted factor P{zi) 
are visibly lower than those with equal-weight. Therefore, 
while approximating the statistical results, the point at Zi is 
endued the weight P{zi) which is correlated to the frequency 
of Zi within the window. 

For the given return value Zi, the frequency N{zi) also 
depends on the length of time windows. Thus, it is worth 
mentioning the influences from the length T on the approx- 
imations. The relations between T and the fitting errors of 
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FIG. 2: (Color online) Log-log plots of weighted fitting errors (ta 
and as vs. T, with T = lOOd ~ 13000d. They exhibit perfect 
power-law behavior The solid lines are power-law regression fits 
over the entire range, which give estimates of power-law exponents, 
rjA = 0.49022 and rys = 0.31122. 




A{z) and B{zY are investigated. To estimate the error purely 
made by window length T, not by information the time series 
embodied, the daily log-returns series is randomly rearranged 
and the fitting errors aA and as vs. T = lOOc? ^ ISOOOd 
are calculated (Fig. |3. It was found that the weighted fitting 
errors a a and ctb of the approximations with weighted fac- 
tor P{zi) decline quickly with a perfect power-law behavior, 
OA = 8.178 X 10^3 . 2^-0.49022 ^jjj ^ ^ 2A% x 10-4 ■ 
y-o. 31122 rpjjjg behavior suggests that the statistical results 
of A{z) (or B{z)'^) with larger T is more feasible to be ap- 
proximated with linear (or parabolic) form. Note incidentally 
that, only the results with the same window length T could 
be compared since different T corresponds to different fitting 
errors. 

The algorithm for the detection of the time-dependence of 
drift term A{z) in this paper can be described as follows: sam- 
ple the long log-returns series with a sliding window of short 
length T, and compute A{z) for each location. The results es- 
timated from the time window which samples a given local 
period, present the corresponding local characters of financial 
markets. The algorithm is more sensitive than merely studying 
transient behavior. The comment for the selection of window 
length T is to choose T long enough so that the averages in Eq. 
Q and Eq. @ are statistically meaningful but not so long as 
to lose the temporal resolution. In the results presented below, 
250(i and 2000d window lengths and bd overlapping (window 
shift by 5 days per time) are used. The corresponding fitting 
errors are, cFA,2bM = 5.46 x 10"'*, (TB,250d = 2.24 x 10"^, 
and (TA,2000d = 1-97 x 10"^, crB,2000d = 1-17 x lO"''. It is 
expected that the variation of A{z) as a function of time can 
accurately indicate interesting dynamical changes in financial 
process. 



III. RESULTS AND DISCUSSIONS 

In the LE the drift parameter A{z) could be seen as an 
action of potential, with A{z) = -W{z). Noted that A{z) 
has a linear form with negative slope, it could be interpreted as 
the effect of linear state-dependent restoring force with sym- 



FIG. 3: The sketch maps of the linear restoring force A{z) and its 
corresponding potential well V{z). 

metrical potential well, 

V{z) = -^z'-bz^f, (8) 

where z presents the position of one given particle enslaved to 
it. The sketch maps of A{z) and V{z) are showed in Fig. |3] 
in which the equilibrium position (A{zq) = 0) is zq — —hja. 
Fig-Elshows the time series of log-returns z(t), restoring force 
A{z) and potential V{z) from 26 May. 1983 to 1 Jan. 1993. 
It is easy to find that the vibrancy of log-returns z presented 
seems to be similar to force A(z) and potential V(z), which 
can be clearly seen from several large events marked in this 
figure, and the restoring force A{z) always presents converse 
effect to log-returns z. In some large events, the potential 
V(z^ is exceedingly large. 

In the following, the time dependence of restoring force 
A^z) will be discussed including the equilibrium position zq 
and the slope coefficient a. In addition, discussions of the dif- 
fusion parameter B(z) and the error analysis will be given. 

A. Equilibrium Position zq 

From Langevin equation JItI . the so called equilibrium po- 
sition Zo, which is the zero value of negative-sloped linear drift 
term, corresponds to the minimum of the potential well. From 
a physical point of view, the average displacement (z) of an 
oscillating particle in the potential well defined by Eq. (jS) 
should also be the minimum of the potential well, (— &/a, 0). 
Thus, we get, 

^0 - (2). (9) 

The average displacement (z) was directly obtained from the 
log-returns series, and the equilibrium position zq was calcu- 
lated from Eq. (|5}. The time dependence of zq and (z) co- 
incide with each other very well all over the ranges [see Fig. 
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FIG. 4: (Color online) (a) Log-returns series z{t), and (b) the corre- 
sponding restoring forces A , and (c) the potential V. Restoring force 
and potential at one given time ti is defined as A{ti) = A{z{ti)) 
and V{ti) = V{z{ti)). The A{z) and V{z) are the average results 
calculated from these T = 250d overlapping windows so far as em- 
bodying the time ti. Some special points are marked with circles to 
reflect relations between the three series. 
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FIG. 5: Sliding window analysis of the time-dependence of (a) the 
average log-return (z), and (b) the equilibrium position zq. The 
lengths of sliding windows are T = 250d (left) and T = 2000d 
(right), and windows shift by 5 days per time, (c) The correlation 
function between zo and {z). Their linear approximations (solid 
lines) are: (z) = 0.9985zo + 3.4060 x 10"" (left) with standard 
deviation SD = 5.78 x 10-^ and (z) = 0.9856zo + 7.8319 x 10"'' 
(right) with standard deviation SD = 1.00 x 10~^. 



|5ja) and|5jb)], and the plots of zq as a function of (z) in Fig. 
ISjc) was excellent agreement with Eq. (|9}. 

In language of finance, the average log-return (z) describes 
the macroscopical trend of the price movement: (z) > in- 
dicates going up, and (z) < indicates going down. There- 
fore, based on Eq. (|9}, the rising trend of prices could be 
estimated if the statistical result zq = —a/6 > 0; otherwise. 



the falling trend could be estimated if zq = —ajh < 0. Thus, 
the equilibrium position zo of restoring force A{z) is compre- 
hended as the 'trend index' of stock prices. Furthermore, the 
stock prices and their log-returns are macroeconomic indica- 
tors which are widely used because of the strong correlation 
between financial markets and economic development. In this 
case, the equilibrium position zq, which is derived from the 
LEs description of financial time series, would be another im- 
portant indicator of macroeconomics. 



B. Slope Coefficient a 

In Fig. |6jc), the time dependence of the slope coeffi- 
cient a calculated from the daily log-returns of S&P 500 with 
T ~ 25Qd and 2000(i are plotted. The ranges of the slope 
coefficient a with T = 250d and 2000d are [-1.127, -0.574] 
and [-1.009,-0.699] respectively, both of which are close 
to the value —1. As shown in Fig. (|3} we know that in a 
certain given position z^, steeper (flatter) slope of A{z) cor- 
responds to larger (smaller) restoring force A{zi). Thus, one 
can imagine the mechanism of our model: it would take few 
times for larger forces (slope: a < — 1) to draw particles from 
one side of the equilibrium position zq to another side, which 
we called 'zQ-crossing' action for the moment, and more times 
for smaller forces (slope: — 1 < a < 0). To discuss the afore- 
mentioned mechanism more accurately, normalized daily log- 
returns are used. 



Z- {z)q 



^/{^')t - {z) 



(10) 



which has zero mean value, {g)T — 0. Here the standard devi- 
ation V of log -returns is defined as the time averaged volatility 
10] and the {■ ■ ■)t denotes an average over the entire length 
of the series within time window T. From Eq. (|9jl, g has the 
equilibrium position zq equaling zero, so that the ZQ-crossing 
action could be reduced to the sign convert of g. Thus the 
mechanism can be described as follows: The sign of g changes 
frequently while the slope is quite steep; on the contrary, same 
signs congregate together and sequences of consecutive '+' or 
'— ' appear when the slope is flat. 

The sign series of daily log-returns has been considered 
to study the daily correlation in log -returns |11]. Those re- 
searchers got the conditional dynamics from the sequence of 
consecutive '+' and ' — '. In this paper, however, we will com- 
pare the sequence of the same sign with the converting sign of 
two neighboring days to check the relationship between a and 
ZQ-crossing. 

The sign-cases of a given day and its previous day are: 

'+- 1-', 'H — ', ' — h' and ' '. One can define '++' and 

' ' as the sign-sustained cases, and 'H — ' and ' — h' as the 

sign-convert cases. Incidentally, the contribution of '++' to 
the sign-sustained cases was shown to be a little more than 

' ' on average over the whole series. Then, the time series 

of signs of g is investigated by counting the frequencies of 
sign-sustained cases, m, and sign-convert cases, n. The time 
dependence of the proportion (3 = m/n and the slope coeffi- 
cient a calculated from the daily log-returns of S&P 500 with 
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slope a larger than —1, and f3 larger than 1, with mean values 
a\T=250d = -0.883, a\T=20Q0d = -0.874, and /3\T=250d = 
1.242, /3|T=20ood = 1.223. These imply that, in practice, 
the flat-slope restoring force is more prevalent than the steep- 
slope one, and the case of sign-sustained is more than that of 
the sign-convert. Consequently, from the general appearance 
of sign-sustained cases, the same conclusion was reached as 
that of 1 1 1]: the return of the price during a given day can be 
correlated with the previous day, in particular with the sign 
of the previous day. However, it is worth noticing that, the 
analysis with Langevin approach is more general because it 
contains the positive daily correlation as the particular case 
with flat-slope restoring force, and non-correlation with steep- 
slope restoring force. 

On the other hand, the autocorrelation function C{At'), a 
typically important statistics of stochastic processes, is always 
used to investigate pairwise correlation of the log-returns of 
a financial asset. In the following, we compared it with the 
slope coefficient a from mathematical relations and statistical 
results. It is known that. 



FIG. 6: (Color online) Sliding window analysis of the time- 
dependence of (a) the proportion of sign-convert cases and sign- 
sustained cases /3, (b) the autocorrelation function CAt=i, and (c) 
the slope coefficient a with At = Id. The lengths of sliding 
windows are T = 250d (left) and T = 2000d (right), and win- 
dows shift by 5 days per time, (d) The correlation function be- 
tween a and C. Their linear approximations (solid lines) respec- 
tively are: C = 0.9942 + 0.9958a (left) with standard deviation 
SD = 3.44 X 10~^ and C = 0.9932 + 0.9938a (right) with stan- 
dard deviation SD = 6.44 x lO""*. 



T = 250d and 2000d are plotted in Figs. Igja) and (c). The 
good similarity of f3 and a proves that the slope coefficient 
is related to the zo-crossing action and reflects the correlation 
of neighboring daily log-returns. In detail, in one sampled 
period, while the slope of the restoring force is flat, the given 
day's sign of log-return is more likely to be the same as its pre- 
vious day; on the country, while the slope is steep, the given 
day's sign is more likely to be different from its previous day. 
Thus, the mechanism of the daily coiTelation in log-returns is 
qualitatively explained by the restoring force. 

It's easy to notice from the time dependence of a and /3 [Fig. 
|6ja)(c)] over the whole series, most of the periods have their 



CiAt') = 



{z{t) ■ z{t + At'))T - {z{t))\ 



(11) 



where (■ • ■)t denotes time averaging over all the trading days 
within the sampled local period with length T, At' is time 
increment. 

For the weighted linear Least-squares fit, which has been 
used to approximate the statistical results of A{z), the 
weighted objective function fT3.fl^ can be written as. 



q^J2W^[y^-ic 



(12) 



where Wi presents the weight of the point {xi ,yi). The values 
of a and b corresponding to the minimal values of function 
q are something to be sought for Thus from Eq. ( I12t . the 
solution of slope coefficient a in y{x) = ax + b is achieved, 

_ E^ W^ Y.^ X^y^W, - J2^ X^W^ J2^ y^W^ 

In the previous discussion, the weight of Zi was defined as its 
probability P{zi) (Eq. 0). Thus one get. 
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(14) 



In this paper, a is calculated by this statistic formula. Substi- tuting Eq. (|3} and Eq. Q into Eq. (I14t . we get, 

^ {z ■ A{z))r - {z)t{A{z))t 
{z')t - {z)l 
1 Az{t) ■ z{t + At))T - {zit))T{z{t + At))T 



At' 



{zitrw - {zm 



(15) 
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When T, compared to At, is sufficiently long, the sta- 
tionary assumption of financial time series is: {z{t))T m 
{z{t + At))rp. Then, compared Eq. il5\ with Eq. il 1> . a 
simple relation between a and C will be found. 



a{At) 



1 

At 



[C(At) - 1] 



(16) 



which is valid for any value of At because of no limit to At 
during the derivation. However, since the lack of correlation 
for At > Id, only the case with At — Id is analyzed. The 
statistical result CAt=id is compared with the slope coefficient 
a, and the correlation function between a and C exhibited the 
good effectiveness of Eq. (I16> [showed in Fig. O. The ana- 
lytical results indicate that the sign-cases, slope coefficient a, 
and autocorrelation function CAt=id reflect the similar prop- 
erties of time series. Known that correlations observed in fi- 
nancial time series show the incompleteness of the efficient 
market hypothesis fTlll . the three coefficients C, a, and (3 may 
probably indicate the degree of market efficiency. From the 
same tendency of the three coefficients showed in Fig. |6lwith 
T = 2000(i (right), two conclusions will be arrived at: (i) the 
market lost efficiency from 1961 to 1976 relatively, since the 
values is much larger than the remaining 25 years; (ii) the mar- 
ket tended to be more and more efficient from 1968 to 1999 
because of the decreasing trend of the value. 

From the preceding analysis of the equilibrium position zq 
and the slope coefficient a, the final form of restoring force 
A{z) will be got by substituting Eq. (|9} and Eq. ( I16> into Eq. 
(0, 



Aiz) ^ —{C - l){z 



(17) 



This new form as a function of the traditional qualities, C and 
(z), is a more direct way to understand the dynamical behavior 
of time series. To the financial data, the information given by 
A{z) mixes features of the macroscopical properties together 
with the detail of the prices evolution: the macroscopical trend 
of prices is presented by the equilibrium position zq, and the 
detail correlation between two neighboring days is exhibited 
by the slope coefficient a. 

Eq. ilH can be informatively rewritten as A{z) = — (z — 
zo)/to. to is the characteristic relaxation time, to — —At/a, 
and reflects the same properties of the log-returns as the slope 
coefficient a does, but to is more visualized. The maximum 
values of to, which were calculated from the sliding windows 
of various length, are all less than 2 days, and the average val- 
ues of to withT = 250dandT = 2000rfareto|25od = 1.153d 
and to|20oo<i — l-157(i. Thus, the effect of the restoring force 
A{z) of one given day decays with the characteristic relax- 
ation time less than 2 days. 



C. The diffusion parameter B{z) 

The diffusion parameter B{z), which usually has the form 
shown in Eq. (|6}, corresponds to a state-dependent linear mul- 
tiplicative noise term B{z)dw in Eq. That is to say, the 
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FIG. 7: (Color online) (a) Daily log-returns of S&P 500 index from 
1950 to 1999. (b) The time-dependence of fitting errors a a and ob 
with T — 250d. (c) The time averaged volatility of the same series 
with T = 250d. 



Langevin description of the log-returns requires a linear mul- 
tiplicative noise term to describe the variability of the log- 
returns, which can be interpreted as the variability of log- 
returns increases with log-returns itself. Thus, we conjecture 
that the heavy tailed probability densities are due to the form 
of diffusion parameter B{z) in our Langevin description. 



D. Fitting error a 

It is valuable to investigate the weighted fitting error of 
A{z) and B{z). Fig. 0b) shows the time dependence of (ta 
and (Tb calculated from the sliding window with T = 250d 
down the log-returns series. The time averaged volatility of 
log-returns, which can measure the degree that the market is 
liable to fluctuate, is calculated from the same 250d sliding 
window [see Fig. Etc)]. Compared with the original daily 
log-returns series [Fig. Eta)], one can easily find that the vari- 
ation of aA and as [Fig. 0^)]^ together with the time av- 
eraged volatility [Fig. Etc)], show sudden jumps when very 
volatile periods enter or leave the time window. For example, 
as pointed out by arrows, the jumps at T = 9500d and 3110d 
are caused by the crashes in May. 1962 and on the 'black 
Monday', 19 Oct. 1987. 

Based on the above analysis, one can conclude that the fit- 
ting error, a a and (Tb, sensitively respond to the volatility of 
the financial markets, i.e., the information of historical events 
in real market would influence the accuracy of the fitting re- 
sults presented in this paper. 
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TABLE I: The results from the log-return series of the S&P 500 
index on the time scale r from Id to 5d. The average slope a, the 
average value of autocorrelation C, the average proportion of the 
sign-sustained and sign-convert cases /3, the characteristic relaxation 
time to with T = 250d and 2000d are listed in this table. 

r(d) a C 13 to{d) 

250d 2000d 250d 2000d 250d 2000d 250d 2000d 

~i -0.883 -0.874 0.115 0.125 1.242 1.223 1.153 1.157 

2 -0.467 -0.462 0.533 0.539 2.273 2.252 2.178 2.187 

3 -0.303 -0.299 0.697 0.702 3.205 3.137 3.395 3.414 

4 -0.224 -0.222 0.776 0.779 4.136 4.059 4.416 4.633 

5 -0.182 -0.177 0.818 0.825 4.986 4.886 5.759 5.822 



E. The results of the log-return series z{t)\r!^id 

In the way of our Langevin approach to the log-return series 
z{t)\T- with T = 2d, 3d, 4d, 5d are also studied. It can be con- 
cluded that, the forms of A(z) and B{z) [Eq. Q and Eq. ©I, 
the correlation function between the equilibrium position zq 
and the average log-return (z) [Eq. the correlation func- 
tion between slope a and autocorrelation C [Eq. (II6I1I. and the 
new form of restoring force A{z) [Eq. (II 7H are all effective 
for these series. Table 1 lists the results from the log-return 
series of the S&P 500 index on the time scale r from Id to 5d. 



IV. SUMMARY AND OUTLOOK 

In this paper, we present a coarse-grain time-dependent 
Langevin description of the dynamics of stock prices, which is 



proved to be effective by the results obtained from analyzing 
the S&P 500 index. The time dependence of drift parame- 
ter A{z), which was considered as the restoring force, was 
investigated by the simple sliding windows algorithm. Signif- 
icantly, while choosing the right weighted factor (Eq. 0) to 
approximate the statistical results of A{z), the hnear approx- 
imations of A{z) can reflect both the macroscopical and the 
detail properties of the price evolution, and the final form of 
the restoring force Eq. (I17> can be achieved from analytical 
methods. The macroscopical trend of price could be investi- 
gated from the equilibrium position zq, and the daily coiTe- 
lation in log-return was exhibited by the flat slope coefficient 
a. The mechanism of our model is discussed by analyzing the 
sign series of log-returns. Therefore, from the restoring force 
A{z) in Langevin approach, one can get the properties of ex- 
perimental data or the properties of financial markets. Further- 
more, it must be pointed out that the random force B{z) also 
plays an important role in the dynamics of financial markets, 
which will be addressed further in the future study. 
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